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ABSTRACT 

Global magnetic field extrapolations are now revealing the huge complexity of the Sun’s 
corona, and in particular the structure of the boundary between open and closed magnetic flux. 
Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs 
in highly fragmented current layers, and that this typically leads to a dramatic increase in the 
topological complexity beyond that of the equilibrium field. In this paper we use static models 
to investigate the consequences of reconnection at the open-closed flux boundary (“interchange 
reconnection”) in a fragmented current layer. We demonstrate that it leads to efficient mixing of 
magnetic flux (and therefore plasma) from open and closed field regions. This corresponds to an 
increase in the length and complexity of the open-closed boundary. Thus, whenever reconnection 
occurs at a null point or separator of this open-closed boundary, the associated separatrix arc of 
the so-called S-web in the high corona becomes not a single line but a band of finite thickness 
within which the open-closed boundary is highly structured. This has significant implications 
for the acceleration of the slow solar wind, for which the interaction of open and closed field 
is thought to be important, and may also explain the coronal origins of certain solar energetic 
particles. The topological structures examined contain magnetic null points, separatrices and 
separators, and include a model for a pseudo-streamer. The potential for understanding both the 
large scale morphology and fine structure observed in flare ribbons associated with coronal nulls 
is also discussed. 

Subject headings: Sun: corona; Sun: magnetic fields; magnetic reconnection; Sun: solar wind 


1. Introduction 

It is well established that the magnetic field 
that permeates the solar corona has a highly com¬ 
plex structure. Although it is very difficult to 
measure directly the magnetic field vector in the 
corona, this complexity can be inferred from ob¬ 
servations of the line-of-sight magnetic field at the 
photosphere. With each new satellite mission that 
is launched, we observe photospheric magnetic 


flux concentrations on ever smaller scales (that 
seem to exhibit a power-law distribution with size, 
Parnell et al .|20 09). Magnetic field extrapolations 


based on these observed photospheric polarity dis¬ 
tributions exhibit an often bewildering degree of 
complexity. Understanding the evolution of such a 
complex magnetic field structure is a major chal¬ 
lenge. 

In recent years, significant progress has been 
made in developing tools with which to charac- 
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terise the coronal magnetic field. One approach 
involves segregating the photospheric magnetic 
field into discrete flux patches. This then allows 
the corona to be divided into distinct domains, 
each defined by the flux connecting pairs of these 
patches. Between these coronal flux domains are 
separatrix surfaces , that emanate from magnetic 
null points. The intersection of two such separa¬ 
trix surfaces forms a separator field line - a field 
line that connects two null points and lies at the in¬ 
tersection of four flux domains. Indeed, magnetic 
field extrapolations reveal the presence of a web 
of null points, separatrix surfaces, and separators 
that form a skeleton based upon which the mag¬ 
netic connectivity of the coronal field may be un¬ 


derstood (e.g. Longcope 2005 Regnier et al. 2008). 
The separatrix surfaces of this skeleton represent 
locations at which the mapping between boundary 
points via the magnetic field lines exhibits discon¬ 
tinuities. Also of interest are layers in which this 
field line mapping exhibits strong (but finite) gra¬ 
dients. These are known as quasi-separatrix layers 
(QSLs), being regions at which the squashing fac¬ 
tor, Q ( |Titov et al.||2002 ), is large. 

Null points, separators, and QSLs, at which the 
field line mapping is either discontinuous or varies 
rapidly, are of interest not only in analysing the 
structure of the coronal magnetic field, but for un¬ 
derstanding its dynamics. This is because these 
locations are prime sites for the formation of cur¬ 
rent layers at which magnetic reconnection may 


occur, releasing stored magnetic energy (Pontin 


2012 and references therein). In particular, they 


have been implicated in the formation of current 
sheets associated with solar flares, jets, and coro- 


nal mass ejections (e.g. Fletcher et al. 2001 

Man- 

drini et al. 

2006; Barnes 2007| Titov et al. 

2008 

Lynch et a 

. 2008 Pariat et al. 2010[ Sun 

et al. 

2013). One major piece of supporting evidence 


is the coincidence of flare ribbons with footpoints 
of separatrix and QSL field lines in coronal field 


2014). 


extrapolations (Masson et al. 2009 Janvier et al. 


One particular location at which the magnetic 
field line mapping is discontinuous is at the in¬ 
terface between closed and open magnetic flux, 
i.e. the boundary between magnetic field lines that 
are anchored at both ends at the photosphere, and 
those that extend out into the heliosphere. Mag¬ 
netic reconnection at this open-closed flux bound¬ 


ary is one of the principal mechanisms proposed 
to explain the properties of the slow, non-steady, 
solar wind (e.g. Fisk et al. 1998). The slow so¬ 
lar wind is characterised by strong fluctuations in 
both velocity and plasma composition, the latter 
of which is consistent with the wind being com¬ 
posed of some component of closed-flux coronal 
plasma (e.g. |Hansteen and Velli 2012). Reconnec¬ 
tion specifically at the open-closed flux boundary 
is also implicated in the generation of impulsive 
solar energetic particle (SEP) events, due to the 


observed ion abundances of these events (Reames 


2013). 


Typically, computational models of the Sun’s 
global magnetic field exclude the outflowing 
plasma of the solar wind, but include its effect 
by imposing a magnetic field that is purely radial 
at some height above the photosphere (termed 
the ‘source surface’). Excluding all contributions 
to the solar magnetic field other than the global 
dipole, the coronal field is characterised by two po¬ 
lar coronal holes of open magnetic field lines and a 
band of closed flux around the equator, these two 
being separated by separatrix surfaces that meet 
the base of the heliospheric current sheet (HCS) 
at the source surface. The question arises: when 
the full complexity of the coronal field is intro¬ 
duced, what is the nature of the boundary between 
open and closed flux? In a series of papers, Fisk 
and co-workers (e.g. |Fisk et al.||1998[ |Fisk and 


Schwadron 2001 Fisk 2005) developed a model 


for the dynamics of the Sun’s open magnetic flux, 
that was also used to explain the acceleration of 
the solar wind mediated by reconnection between 
open and closed field lines (termed interchange 
reconnection by Crooker et al. [2002). In their 
model, open field lines can freely mix with and 
diffuse through the closed field regions, and in¬ 
deed it is predicted that this open flux component 
should become uniformly distributed throughout 


the (predominantly) closed field region (Fisk and 
Zurbuchen |2006). While noting that such a sce¬ 


nario requires the presence of current sheets in the 
corona between open and closed flux, these stud¬ 
ies do not address the magnetic field structure 
in detail. Indeed, the topological admissibility of 
such free mixing of open and closed flux has since 


been questioned (Antiochos et al. 2007), making it 


difficult to reconcile the interchange reconnection 
solar wind acceleration mechanism with the broad 
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observed latitudinal extension of the slow solar 
wind streams (up to 60°, especially at solar min¬ 


imum) (e.g. McComas et al. 20001. Nonetheless, 
recent modeling of the global coronal magnetic 
field has suggested a resolution to the apparent 
contradiction that plasma that appears to orig¬ 
inate in the closed corona is observed far from 
the HCS at large radii. It has been demonstrated 
that additional regions of open flux that are dis¬ 
connected from the polar coronal holes (at the 
photosphere) may indeed exist. The distinct pho- 
tospheric regions of open magnetic flux are par¬ 
titioned by multi-separatrix structures associated 
with multiple nulls points, typically comprising a 
dome-shaped separatrix enclosing the closed flux 
between the two open field regions, intersecting 


with a vertical separatrix curtain (Titov et al. 


2011 Platten et al. 2014). Even when coronal 


holes are not disconnected, there may exist very 
narrow channels of open magnetic flux at the pho¬ 
tosphere connecting two larger open flux regions. 
In this case the narrow channel is associated with 
a QSL curtain. Both the QSL and separatrix cur¬ 
tains extend out into the heliosphere, and have 
been shown to map out a broad latitudinal band 
around the HCS, termed the S-web (Antiochos 
et al.pOlT Crooker et al.|[2012 ). The correspond¬ 


ing arc structures at the source surface in global 
models are associated with pseudo-streamers in 
the observations, and there is growing evidence 
that these structures are associated with slow so¬ 
lar wind outflow (e.g. Owens et al. [20131. 


The above studies have revealed that the open- 
closed flux boundary has a complex topological 
structure involving null points and their associ¬ 
ated separatrices, separators and QSLs. Moreover, 
it has been recently demonstrated that when re¬ 
connection occurs in astrophysical plasmas, the 



Fig. 1.- A model magnetic field showing the mag¬ 
netic topology following the initial bifurcation of 
the central 3D null following the tearing instability 


- after Wyper and Pontin (2014a). 


the topological effect of reconnection in configura¬ 
tions defined by an isolated separatrix dome and 
a separatrix curtain. In Sections [5] and |6j respec¬ 
tively, we present a discussion of the results and 
our conclusions. 

2. Non-linear tearing during 3D reconnec¬ 
tion 

In recent years a major advance in our under¬ 
standing of magnetic reconnection has been the 
realisation that the reconnection rate can be sub¬ 
stantially enhanced when the current layer breaks 
up in response to a tearing instability. While the 
linear phase of the classical tearing mode is slow, 
non-linear tearing in two dimensions (2D) via the 
plasmoid instability can grow explosively, and lead 
to a reconnection rate that is only weakly depen- 


3D topological complexity can dramatically in- 

dent on the resistivity I 

Loureiro et al. 2007 

Bhat- 

crease beyond that of the equilibrium field (e.g. 

tacharjee et al. 2009 

Loureiro et al.||2013) 

In 2D 


In this paper we use simple static magnetic field 
models to investigate the implications of recon¬ 
nection for the magnetic field connectivity at the 
open-closed flux boundary when the reconnecting 
current layer exhibits a fragmented structure, ex¬ 
pected to be the typical case in the corona. The 
paper is structured as follows. In Section [2] we 
summarise recent relevant results on current layer 
instabilities. In Sections [3] and [3] we investigate 


islands are formed as current layers fragment into 
chains of X-type and O-type nulls. The condi¬ 
tions for onset of the instability are that the (in¬ 
flow) Lundquist number be > 10 4 , and the current 
layer aspect ratio be > 50. 

It has recently been demonstrated that the plas¬ 
moid instability also occurs in three-dimensional 


(3D) current layers. Daughton et al. (2011) per¬ 


formed 3D particle-in-cell simulations of an ini- 
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tially planar, infinite current layer. They noted 
that magnetic ‘flux ropes’ formed in place of the 
magnetic islands from the 2D picture with flux of- 

multiple flux ropes, 
by contrast studied 
MHD simulations of a 3D magnetic null point un¬ 
dergoing external shear driving. They observed 
the initial formation of a laminar current layer 
centred on the null, which was found to become 
unstable at a threshold similar to the 2D case 
(Lundquist number > 2x 10 4 , aspect ratio > 100). 
This threshold onset condition is very likely to be 
exceeded for typical current sheets formed in the 
corona. In this study, it was demonstrated that 
the onset of tearing leads to the creation of new 
3D null points in bifurcation processes. In par¬ 
ticular, 3D spiral nulls are formed that are the 
analogue of 2D islands, and the spine lines of each 
of these nulls forms the axis of a pair of magnetic 
flux ropes, as shown in Figure [l] Crucially, and in 
contrast to the 2D case, these flux ropes are open 
structures - they are not surrounded by flux sur¬ 
faces as 2D islands are (this is required on topolog¬ 
ical grounds due to the variation of B along the di¬ 
rection of the flux rope axis and the solenoidal con¬ 
dition on B). Asa result, no new isolated domains 
of magnetic connectivity are formed. Rather, the 
flux ropes are composed of a mixture of flux from 
the two connectivity domains (flux located above 
and below the separatrix of the single null prior to 
the instability), wrapped around each other. The 
result is that, when the field line connectivity is 
analysed, a layer is found in which magnetic flux 
from the two connectivity domains is rather effi¬ 
ciently mixed in spiral patterns associated with 
the multiple flux rope structures. Our intention 
here is to investigate the effect on the global con¬ 
nectivity when the fragmented, reconnecting cur¬ 
rent layer is embedded in some generic coronal 
field structures. 


ten threading in and out of 
Wyper and Pontin (2014a|b 


3. Isolated dome topology 


We first examine the simplest generic coronal 
configuration containing a 3D null point the iso¬ 
lated dome topology (see Figure [2](a)). Such a 
configuration is always present, for example, when 
a photospheric region of one magnetic polarity is 
embedded in a region of opposite polarity. Such 
3D nulls are preferred sites of current sheet forma¬ 


tion (Pontin and Craig 2005 Pontin et al. 2007), 




Fig. 2.— (a) Magnetic field lines for the isolated 
dome topology given by Equation ([!]), and (b) af¬ 
ter the addition of one flux ring (state la) - sam¬ 
ple field lines in the flux ropes are coloured red. 
Green boxes outline the regions in which connec¬ 
tivity and Q maps are calculated in Figures [3j |4j 
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and the significance of reconnection in an isolated 
dome configuration has been considered, for exam¬ 


ple, by Antiochos (1996); Pariat et al. (20091; Ed¬ 


mondson et al. (2010). As demonstrated by Pon- 


tin et al. (2013), spine-fan reconnection in such 


a dome topology is characterised by a transfer of 
flux in one side of the separatrix dome and out the 
other side, which can be driven dynamically or oc¬ 
cur during a relaxation process as the coronal field 
seeks a minimum energy state. 

In order to examine the effect of tearing of 
the null point current sheet on the field struc¬ 
ture we consider the following simple model. The 
fields that we construct are not equilibrium fields 
(e.g. are not force-free) - however this is of no 
importance since we do not consider here any dy¬ 
namical processes. Rather, our purpose is solely to 
examine the field topology/geometry that results 
from breakup of a reconnecting current layer. We 
consider the (dimensionless) field 

B = xe x +ye y -|-(l — 2z — 4r 2 tanh 2 (4r)) e z , (1) 

where r = \Jx 2 + y 2 . This magnetic field contains 
a null point at ( x,y,z ) = (0,0,0.5) above a pho- 
tospheric plane represented by z = 0, see Figure 
j2ja). The isolated dome topology appears over a 
wide range of scales in the corona. Hereafter we 
discuss any length scales in terms of a characteris¬ 
tic ‘macroscopic’ length scale of the overall struc¬ 
ture, that we denote D. In our model field both 
the separatrix dome diameter at z = 0 and the 
null point height are of order 1, so for the model 
field D ~ 1. From magnetic field extrapolations 
it is observed that in the corona this scale D of 
the dome separatrix structure can be as large as a 
few hundred Mm (usually in the vicinity of active 
regions, e.g. Del Zanna et al.||2011 Platten et al. 


2014), and at least as small as a few tens of km 


(in quiet-Sun regions, where the lowest null point 
height in extraploations is likely limited by the 
magnetogram resolution, e.g. Regnier et al.|20 08). 

Onto the ‘background’ field |l]) we super¬ 
impose a magnetic flux ring to simulate the topo¬ 
logical effect of magnetic reconnection occurring 
at some particular location in the volume. This 
method is motivated and described in detail in 


Wilmot-Smith and Hornig ( 

2011 

); 

Pontin et al. 


(2013). The field of this flux ring is taken of the 


form 

B r = B 0 V x (A y e y ), 

(■ x-x N ) 2 ( y-yN ) 2 (z-z N ) 2 


A y = exp - 


L 2 


l 2 


L 2 


We begin by super-imposing a single flux ring of 
this form onto the field of Equation (jT|) , centred 
at the null (i.e. we set xn = yN = 0, zn = 0.5). 
L and l are the characteristic size of the flux ring 
in the xz -plane (plane of B#) and along y, re¬ 
spectively. For Bq small, the effect of adding the 
flux ring is to collapse the spine and fan of the 
null point towards one another, as described by 


Pontin et al. (2013). This has the effect of trans¬ 


ferring flux in one side of the dome and out of the 
other, and is consistent with the topology of a sin¬ 


gle laminar reconnection layer at the null (Pontin 
et a!. 12007) Galsgaard and Pontin|20TT ). However, 


for larger values of B 0 the field becomes elliptic 
at ( x,y,z) = (xN,yN, zn) as the flux ring field 
dominates over the hyperbolic background field. 
The result is a bifurcation of the original null into 
three null points as described above in Section [2] 
see Figure [l] - and the generation of a pair of flux 
ropes (see Figure[2](b)). This models the magnetic 
topology when the current layer undergoes a spon¬ 
taneous tearing instability as observed by |Wyper| 


and Pontin (2014a). Parameters for the magnetic 


field - denoted state lb - are presented in Table 

□ 

In order to visualise the new field structure cre¬ 
ated after tearing onset we first plot a connectiv¬ 
ity map of field lines from the lower boundary. 
That is, we trace field lines from a grid of foot- 
points on the lower boundary and distinguish field 
lines that are closed (return to the lower boundary 
z = 0) and open (exit through the upper boundary 
z > 1). The resulting map can be seen in Figure 
[3j While the flux ropes are approximately circular 
near the apex of the dome, they are compressed 
towards the separatrix and stretched in the az¬ 
imuthal direction by the global field geometry, and 
thus appear as flattened spiral structures in the 
connectivity map (Figure [3]) a)). In order to more 
clearly visualise this structure we reproduce the 
map in a polar coordinate system in Figure [3](b) . 
The observed spiral pattern of mixing of open and 
closed flux reproduces the behaviour in the dy¬ 
namic MHD simulations - see Figure 8 of |Wyper| 
and Pontin (2014aI. One can include the effect of 


5 

































































(b) 


0.730 

0.728 

0.726 

0.724 

0.722 

0.720 

0.718 



0.730 
0.728 
0.726 
0.724 
0.722 
0.720 
0.718 

0.716 

0.0 0.5 1.0 

(c) 


1.5 2.0 2.5 3.0 

theta 



Fig. 3.— Pliotospheric connectivity maps for 
the isolated dome topology. White regions con¬ 
tain footpoints of open field lines, black of closed 
field lines (grey: field lines from these points not 
traced), (a) State lb, and (b) the same plot in 
standard r — 8 coordinates centred at (x, y) = 
(—0.01, 0). (c) State lc containing three flux rings. 


a further breakup of the current layer through the 
inclusion of additional flux rings. Adding two such 
flux rings centred at the non-spiral nulls of state la 
leads each of these nulls to undergo a bifurcation, 
resulting in a total of seven nulls. This naturally 
introduces additional spiral structures in the field 
line mapping, as show in Figure [3](c) (state lc, 
for parameters see Table [lj , and if one were to 
iterate this procedure by adding more flux ropes 
a mapping with complexity of the order of that 
seen in the MHD simulations of|Wyper and Pom 


tin (2014a I could be obtained. 


We now turn to consider the characteristics of 
the open flux that exits the domain through the 
top boundary at z > 1. Since all of this flux 
is open, a connectivity map does not reveal this 
structure. However, we can use for example the 
squashing factor , Q , to visualise the field line map¬ 
ping from z = 0 to z = 1. Here we plot Q on the 
surface z = 1. The distribution of Q is obtained 
by integrating field lines from a rectangular grid of 
typically around 10 b footpoints and then calculat¬ 
ing the required derivatives using finite differences 
over this grid. Q is formally infinite on spine and 
fan field lines since they represent discontinuities 
in the field line mapping. However, calculating 
Q numerically as we do here they show up only 
as sharp points and lines, respectively, with very 
high values of Q. One should therefore not at¬ 
tach physical meaning to the maximum value of 
Q in the plots (attained at the separatrix/spine 
footpoints) as it is determined entirely by the res¬ 
olution of the field line grid. 

In the background dome topology of Equation 
0 a single spine line intersects the z = 1 bound¬ 
ary, and the Q-map displays a single maximum 
at the origin. When a null point bifurcation oc¬ 
curs during reconnection, the topological structure 
changes to that shown in Figure [f] In this case a 
vertical separatrix extends up to the top bound¬ 
ary, bounded on either side by a pair of spine lines. 
Examining the Q distribution on z = 1 for state 
la (Figure |4^,), the separatrix footprint is clearly 
in evidence (horizontal line of high Q). Increasing 
the strength of the flux ring in the model (state 
lb) leads to a lengthening of this separatrix due to 
the increased separation of the nulls (Figure |4ju) . 
A further breakup of the current sheet leads to 
the appearance of multiple vertical separatrices, 
as seen in Figure Qc). 
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Fig. 5.— (a,b) Q on the top boundary of the simulation of Wyper and Pontin (2014b) at two times soon 
after tearing occurs, (a) t = 13.8. (b) t = 14.4. (c) Magnetic field line structure at t = 14.4. Blue field lines 
show one pair of flux ropes. The grey iso-surface shows the cylindrical surface used to calculate Q. Volume 
shading indicates the current density modulus. 


There are two additional noteworthy features 
of the Q distributions. First, note the arcs of high 
and low Q that run parallel to the separatrix foot¬ 
print. These become more pronounced and numer¬ 
ous as the flux ring strength is increased (compare 
Figures |4](a,b)). Their origin can be understood 
as follows. Consider held lines traced down from 
the top boundary that enter one of the flux ropes. 
Some local bundles of held lines will spiral around 
the rope axis and then ‘leave’ the flux rope at its 
top or bottom (in z) with a range of values of 6 but 
roughly constant z values (0 being the azimuthal 
angle in the xy- plane). Due to their range of 9 
values on leaving the rope they diverge in the az¬ 
imuthal direction as they are traced onwards to 
the lower boundary - and therefore exhibit rela¬ 
tively high Q values. By contrast, adjacent held 
lines that leave the hux rope along its sides at ap¬ 
proximately equal 9 values but differing z values 
are naturally squeezed in towards the fan as they 
approach the photosphere - they do not diverge 
with the null point fan geometry owing to their 
close alignment in the 9 direction. This leads to a 
lower Q. For a stronger hux ring (more substan¬ 
tial hux rope) held lines have the opportunity to 
spiral multiple times around the rope axis, leading 
to multiple Q stripes. The second feature to note 
are the high-Q ridges emanating from each spine 
footpoint in the y-direction, that are present for 
the following reason. Adding the hux rings nat¬ 
urally generates a strong held component in the 


y = 0 plane. Therefore the two non-spiral nulls 
have quite asymmetric fan eigenvalues (their ratio 
is around 2.5 in state la). The weak held direc¬ 
tion corresponds to the y-direction, and it is nat¬ 
ural that Q is largest in this weak-held region of 
diverging fan held lines. 


To test the robustness of the structures in Q 
described above, two Q-maps were calculated us¬ 
ing magnetic helds taken from the dynamic MHD 
simulation of Wyper and Pontin (2014b). To 
avoid discontinuities in Q brought on by the cor¬ 
ners of the domain we calculated Q using the foot 
points of held lines traced from a hxed grid on the 
top boundary to a cylindrical surface dehned by 
r = \Jy 2 + z 2 = 3.4, red and grey surfaces in Fig. 
[5](c) respectively. Fig. [5])a) shows Q on the top 
“open” boundary soon after tearing has occurred 
in the layer - see also Fig. 5 of Wyper and Pon- 


tin 


(2014b). Note that x is the vertical direction 
in the simulation domain and that the spine-fan 
collapse occurs in the z = 0 plane, Fig HH- At 
this time there is a single pair of hux ropes within 
the current layer. Despite significant fine struc¬ 
ture (likely resulting from turbulent dynamics in 
the outflow region) the structures in our simple 
model described above are clearly evident also in 
the dynamic MHD simulation. A short high-Q line 
corresponding to a vertical separatrix surface is 
apparent near (y,z) ~ (—0.935,0), whilst a num¬ 
ber of parallel stripes of Q can be seen extending 
to either side of it. Additionally, two high-Q ridges 
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Fig. 4.— Plot of Q on the top boundary 2=1, 
for (a) state la, (b) state lb, and (c) state lc. 


emanate from the ends of this separatrix surface. 
The observed closed loop of high Q results from 
the flux rope pair being located in the reconnec¬ 
tion outflow having detached from the open-closed 
boundary, see the discussion of|Wyper and Pon- 


tin 


( 2014a|b ). This splits the field lines that con¬ 
nect from x = 0.5 to y = —3.5 into two bundles: 
those that connect directly from the top bound¬ 
ary to the side and those that loop first around 
the back of the flux rope pair. The foot points 
of the latter are found within the loop of high Q. 
Note that the field line connectivity changes con¬ 
tinuously around the boundary of this loop, so the 
value of Q is large but finite. At the later time two 
pairs of flux ropes are present in the outflow re¬ 
gion of the current layer resulting in an additional 
separatrix footprint being present, Fig. [5^b). The 
gap observed between the pair of vertical separa¬ 
trix footprints (in contrast to Fig. |4|c)) is again 
a result of the detachment of the nulls from the 
open-closed boundary. 

We conclude that tearing of the reconnecting 
current layer at an isolated coronal null separa¬ 
trix dome leads to the formation of an envelope 
around the initial dome structure in which mag¬ 
netic flux from inside and outside the dome is effi¬ 
ciently mixed together. Additionally, vertical sep¬ 
aratrix curtains are formed during each null point 
bifurcation. The implications of these results will 
be discussed in Section [5] 


4. Separatrix curtain / pseudo-streamer 
topology 

4.1. Magnetic field model 

Isolated separatrix dome structures associated 
with a single null as considered in the previous 
section separate small pockets of closed flux from 
the open flux in the polar regions (as well as being 
prevalent in closed flux regions). However, it is 
also typical to have much more complicated sep- 
aratix configurations separating open and closed 
flux. In particular, in global field extrapolations 
it is seen that vertical separatrix curtains lie be¬ 
tween coronal holes that are of the same polarity 


but are disconnected at the photosphere (Titov 


et al. 2011 Platten et al. 2014). These curtains 


together with QSLs associated with narrow corri¬ 
dors of open flux, are associated with arc struc¬ 
tures at the source surface in global models that 




















are interpreted as being associated with pseudo 
streamers (Antiochos et al.|2011 Titov et al. 2011 


Crooker et al. 12012). 


We consider here a simple model containing 
a vertical separatrix surface representing one of 
these curtains. This intersects a separatrix dome 
associated with three coronal null points along two 
separator lines. The separatrix curtain consists of 
the fan surface of one of these nulls (see Figure [6]). 
The magnetic field expression for our model is as 
follows 


B = -2xe x - y{y 2 - l)e y 

+ {(z — 1)(1 + 3 y 2 ) + 2x 2 
—2(tanh(6y + 8.4) — tanh(6j/ — 8.4) — 2)y 2 } e,. 

(3) 


Again a characteristic length scale of the overall 
structure is of order 1 in the model field (say the 
null point height or separation - see Figure [b]), 
that we refer to as D. On the Sun, D is ob¬ 
served to be as large as a quarter of the solar 
radius (~ 170Mm, e.g. Wang et al. 2007), and 


may be at least as small as tens of kilometers in 
quiet Sun regions as discussed before. Note that to 
find the separators in these models the “progres¬ 
sive interpolation” method (Haynes and Parnell 


2010 Close et al. 2005) was used, whereby field 


lines were traced from a ring encircling one of the 
associated null points on its fan plane to identify 
the approximate position of each separator, before 
using an iterative bisection procedure to find each 
separator to a desired accuracy. 

4.2. Separator current layer breakup 

4-2.1. Magnetic field connectivity 

Connecting the three null points along the top 
of the separatrix dome are a pair of separator field 
lines, Fig. [6| Like 3D null points, these are known 
to be preferred sites for current sheet formation 


and magnetic reconnection (e.g. Longcope 1996). 
It has been previously observed that these sheets 
are prone to fragment, yielding a current layer con¬ 
taining multiple separators (|Parnell et al. 2008 


2010). We thus begin by considering the effect of 


super-imposing one then more flux rings to sim¬ 
ulate the effect of tearing in a reconnecting cur¬ 
rent layer around the separator. Such a model 
with a single flux ring was presented by |Wilmot- 


Smith and Hornig (2011) using a background field 


with two nulls, both with initially planar fan sur¬ 
faces. They noted that as they increased the 
strength of the flux ring new separators appeared, 
coinciding with the formation of distinct new do¬ 
mains of magnetic flux connectivity. By magnetic 
flux domain here and throughout we mean a vol¬ 
ume within which there is a continuous change of 
field line connectivity. Distinct flux domains are 
bounded by separatrix surfaces at which this con¬ 
nectivity change is discontinuous. 


We observe the same effect when adding a single 
flux ring on the separator - state 2a (see Table [TJ. 
Tracing field lines from the photosphere (z = 0) 
and making a connectivity map as before, we ob¬ 
serve the presence of a region of open flux nested 
within the closed field region, Figure [7](a,b). In¬ 
creasing the strength of the flux ring, we observe 
progressively more open and closed flux volumes 
being created, nested within one another - Fig. 
[7](c) , state 2b. The formation of these new flux 
domains corresponds to the formation of new pairs 
of separators joining the two associated nulls. Fig¬ 
ure |8](a) demonstrates this for state 2b. Whereas 
originally one separator joined the central and end 
nulls the formation of the three nested flux do¬ 
mains (Fig. [7](b)) corresponds to the birth of three 
additional pairs of separators, giving seven in to¬ 
tal (green and purple field lines), see below for a 
further discussion. 


We now turn to examine the connectivity of 
field lines that extend outwards into the helio¬ 
sphere (those that exit through the top boundary). 
Throughout we consider the surface z = 2 as being 
the ‘top’ boundary - field lines are close to vertical 
above this plane and so little deformation of the 
field line mapping occurs. In Figure [9])a) a map 
of Q is plotted on the top boundary (as calculated 
between the two surfaces z = 0 and z = 2) for 
state 2b. Note that a colour scale is not shown 
since the maximum value is arbitrary, depending 
only on the resolution of the field line grid. We 
observe the imprint of the separatrix curtain, as 
well as additional nested loop structures that cor¬ 
respond to additional separatrix surfaces separat¬ 
ing nested flux domains. In Figure [9](b) a connec¬ 
tivity map is plotted - field lines that intersect the 
black region connect to the photosphere at x < 0 
on one side of the separatrix dome/curtain struc¬ 
ture, while field lines intersecting the white region 


9 





































Fig. 7. - Photospheric connectivity maps for the separatrix curtain topology. White regions contain foot- 
points of open field lines, black of closed field lines, (a) State 2a, (b) state 2a, close-up, (c) state 2b. 



Fig. 6.— Magnetic field lines for the separatrix 
curtain topology, Eq. ([3|. Red and blue spheres 
correspond to nulls with topological degree of +1 
and —1 respectively, whilst the two separators are 
shown in green and purple. Green boxes outline 
the regions in which connectivity maps are calcu¬ 
lated. 


connect to the photosphere on the other side of the 
dome, x > 0. Embedding our structure in a global 
field these two different regions would correspond 
to open field regions of the same polarity that are 
disconnected at the photosphere (see e.g. Figure 5 
of Platten et al.||2014 ), and thus the figure shows 
that flux from the two disconnected coronal holes 
forms a mixed, nested pattern. These nested con¬ 
nectivity regions are entirely equivalent to those 
described in the photospheric connectivity maps 
above. It is expected that in a dynamic evolution 
there is continual reconnection of field lines within 
the current layer, and thus a mixing of plasma be¬ 


tween all of the nested flux domains (see e.g. Par- 
nell et al.~p 0081. As such, field lines at large height 


in these nested flux regions will continually be re¬ 
connected with those from the closed flux region. 


The addition of further flux rings - represent¬ 
ing further plasmoid structures in the reconnect¬ 
ing current layer - leads to the formation of ad¬ 
ditional adjacent sets of nested open/closed flux 
domains. Figure [9](c) shows the connectivity map 
when three flux rings are present, state 2c. The 
complexity quickly becomes very high, with ex¬ 
tremely thin layers of connectivity with character¬ 
istic thickness of order 10 -3 .D - even though the 
flux rope structures and their collective footprint 
in the solar wind remain much larger, having di¬ 
ameters of order 10 _1 D. The inclusion of further 
flux ropes would decrease the length scales in the 
mapping yet further. This complexity of field line 
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Fig. 8.— (a) State 2b. Middle panel shows the 3D structure of the magnetic field. Red and blue spheres 
correspond to nulls with topological degree of +1 and —1 respectively. The separators are alternately shown 
in green and purple representing their pairwise formation. Left panel: null graph showing the connectivities 
of the nulls and separators. Right panel: connectivity map taken from a cut at y = —0.5. Yellow and white 
regions correspond to open held lines that connect to x < 0 and x > 0 on the photosphere, respectively. 
Similarly, grey and black regions correspond to closed held lines that connect from the vicinity of the spine 
footpoint at ( x,y,z ) = (0, —1,0) to the vicinity of the spine footpoints of the central null at x < 0 (grey) 
and x > 0 (black). The separatrix surfaces of the two nulls (red and cyan lines) lie at the intersections of 
the coloured regions. A, B and C denote the flux domains corresponding to the flux tubes shown in Figure 
10 The green and purple dots denote where the separators cross this plane, (b) State 3a. Null graph and 
3D held structure when an end null is bifurcated, (c) State 4. Null graph and 3D held structure when 
the central null is bifurcated. Note that the connectivities shown in each null graph neglect the winding of 
different separators around one another. 
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Fig. 9.— (a) Squashing factor Q on the surface 
2 = 2, and (b) close-up connectivity map, both 
for state 2b. The black region contains footpoints 
of field lines that connect to the photosphere at 
x < 0 (on one side of the separatrix dome/curtain 
structure), in the white regions field lines connect 
to x > 0 (on the other side). The letters corre¬ 
spond to the letters marking the flux tubes in the 
3D plot of Figure 10. (c) Connectivity map for 
state 2c. 


mapping was observed in a related context by|Par- 


nell et al. (2010). 


4-2.2. Relation to 3D magnetic topology 


The direct association between the newly- 
created nested flux domains and additional sep¬ 
arators that form in the domain is demonstrated 
in the right-hand frame of Figure [8](a). Here we 
note that the seven separators lie at the inter¬ 
sections of the four different connectivity regions. 
The nested formation of flux regions and the as¬ 
sociated pairs of separators may be understood as 
follows. Consider a bundle of field lines passing 
in along the open spine of the null at (0, —1,1). 
As the flux ring strength is increased some of 
these field lines are wrapped repeatedly around 
the axis of the flux rope before they reach the 
photosphere. Field lines reach the photosphere 
near spine footpoints of the central null, either at 
x < 0 or x > 0. They may do this directly, or by 
first winding once, twice, or more times around 
the flux rope axis. This is demonstrated in Fig- 
where flux tubes are plotted from each 


ure 
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of the nested connectivity domains that intersect 
the upper boundary of state 2b (marked ‘A’, ‘B’ 
and ‘C’ in Figure |)](b)). Each additional winding 
corresponds to a new flux domain. This is be¬ 
cause as the strength of the flux ring is increased 
the separatrix surfaces of the two nulls (red and 
cyan curves in the right panel of Figure |8](a)) fold 
over to intersect with one another multiple times. 
Each additional pair of intersections correspond to 
a pair of new flux domains (bounded by portions 
of the separatrices) and a pair of new separators 
(see the right-hand image of Figure [8|(a) a nd Fig¬ 
ures 3 and 4 of Wilmot-Smith and Hornig 2011). 


Field lines wind more times closer to the flux ring 
axis, and so new topological domains are formed 
within the previous ones along with a pair of sep¬ 
arators. This explains the nested nature of the 
connectivity domains observed in Fig. [7] and why 
each new pair of separators have one half twist 
more than the preceding two, Fig. [81(a). 


4.3. Effect of null point bifurcations 

So far we modelled the case in which the current 
layer forms along the separator, with this current 
layer breaking up but the number of nulls remain¬ 
ing fixed - that is, no null bifurcations occurred. 
However, another distinct possibility is that the 
current layer that forms in response to a dynamic 
driving of the system contains one or more of the 
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Fig. 11.— (a) State 3a: Close-up connectivity map on the photosphere with colours as in Figure 0 (b) 
State 3a: Q-map on the top boundary, (c) State 3a: Close-up connectivity map on the upper boundary with 
colours as in Figure [9j (d) Same as (a) for state 3b. (e) Same as (c) for state 3b. 


coronal nulls. When such a current layer breaks up 
we would expect a bifurcation of the corresponding 
null point (s), which naturally should also coincide 
with the formation of additional separators - in¬ 
deed this may well occur even when the current 
is focussed away from the nulls, as observed by 


Parnell et al. (2010). 


We consider here two distinct cases, in the first 
of which the null point initially at (x, y, z) ss 
(0, —1,1,) is bifurcated into multiple null points, 
and in the second of which the central null (ini¬ 
tially at (0,0,1)) is bifurcated. Consider first the 
situation where the end null is bifurcated. First, 
adding a single flux ring of sufficient strength at 
the initial location of the null we obtain a bifur¬ 
cation to form three nulls as in Figure [I] (state 
3a) - see Figure |8](b). Let us now examine the 
result for the magnetic flux connectivity, consid¬ 


ering first flux intersecting the photosphere. In 
Figure [7T)( a) we see that the photospheric connec¬ 
tivity map appears as before: new nested open 
and closed flux domains are created in the vicin¬ 
ity of the spine footpoints of the central null. The 
connectivity map for open flux traced from the 
upper boundary is shown in Figure ©c) (where 
the colours have the same meaning as before). 
As shown in the Q-map (Figure & , the main 
separatrix curtain is diverted in the positive re¬ 
direction for negative y. It terminates on the spine 
of one of the null points located in the vicinity of 
(, x,y,z) = (0,—1,1). There is then an additional 
separatrix footprint orthogonal to this bounded by 
the spines of the newly created nulls as in Figures 
[T] [fj Interestingly though, the arcs of high Q em¬ 
anating from this separatrix (as in Figure [4]) now 
form the boundaries of the flux domains that con¬ 
nect to opposite sides of the dome footprint. This 


13 




















Fig. 10.— Flux tubes (red, green, magenta) lo¬ 
cated in the nested flux domains that intersect 
the upper boundary of state 2b. The letters corre¬ 
spond to the letters marking the flux domains in 
Figure |9])b). 


can be understood by considering that each arc of 
high Q in Fig. [4] represents a further half turn of 
field lines along the axis of the flux rope, i.e. the 
outermost two arcs correspond to field lines ex¬ 
iting along one or other spine of the central null 
having wound up to once around the flux rope 
axis, the next pair to field lines that first wind 
between once and twice around the flux rope axis, 
and so on. When these field lines are mapped on to 
the photosphere as in the domed single null case 
this leads to the continuous but rapid change in 
connectivity denoted by the Q ridges. When such 
field lines separate along the spine of a distant null 
the change in connectivity becomes discontinuous, 
forming the nested flux domains, see also Section 
l4~2T2l 


As before, the connectivity maps quickly be¬ 
come significantly more complex when additional 
flux ropes are added. Figures 11 'd,e) show the 
photospheric and upper boundary connectivity 
maps when an additional two flux rings are added 
to generate a bifurcation to a state with three flux 
ropes pairs and seven nulls, state 3b. Again, char¬ 
acteristic length scales of the mapping layers of or¬ 
der 10~ 3 .D or below are observed within a mixed 



_ 03 0 3 0.00 0.04 0.08 

(b) x ( c ) 

Fig. 12.— All for state 4: (a) Close-up connec¬ 
tivity map on the photosphere with colours as in 
Figure[7j (b) Q on the upper boundary, (c) Close- 
up connectivity map on the upper boundary with 
colours as in Figure [9] note that the ^-direction 
is stretched in (c) for clarity. 


flux region with dimensions of order 10 1 D. 

Finally, suppose that the fragmentation of the 
coronal current layer leads to a bifurcation of the 
central null point. Adding a single flux ring leads 
to a bifurcation to a state with three null points 
(state 4) as before. Analysis of the resulting topol¬ 
ogy reveals a situation that mirrors state 3a. As 
shown in Figure [l2| a), this time the photospheric 
connectivity maps show adjacent crescent-shaped 
domains of open and closed flux, symmetric about 
y = 1 since the null bifurcation now leads to a 
bifurcation of both of the initial separators. Cor¬ 
respondingly, nested flux domains of alternating 
connectivity are now observed in the connectivity 
map for the upper boundary (Figure |T2|)b) ), this 
time emanating from the footpoints of both of the 
vertical open spines. 
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5. Discussion 


5.1. Mixing of open and closed flux 


As shown in the above models, reconnection at 
the Sun’s open-closed flux boundary can result in 
that boundary taking on a highly non-trivial struc¬ 
ture. In the presence of an isolated null point sep- 
aratrix dome no new flux domains are created but 
an envelope forms around the initial dome struc¬ 
ture in which magnetic flux from inside and out¬ 
side the dome is efficiently mixed together in spiral 
patterns. As shown by Wyper and Pontin (2014a) 
magnetic flux is continually and recursively recon¬ 
nected from open to closed and back again within 
this envelope (i.e. is reconnected back and forth 
multiple times between open and closed regions - 
The result for the field at 


Parnell et al. 

2008 

large heights 

s that 


the original spine line within which field lines are 
being continually reconnected with those from the 
closed region beneath the dome. 


If we consider a more complicated structure in 
which coronal separators are present, the breakup 
of the current layer leads to the formation of new 
flux domains. In particular, new open and closed 
magnetic flux domains form in nested structures, 
whose length scales become rapidly shorter; even 
for the models considered here containing just 
three flux rope pairs characteristic length scales 
of the mapping layers of order 10~ 3 D or smaller 
are observed. The expectation is that in a dy¬ 
namic evolution, continual transfer of flux/plasma 
between the narrow open and closed layers would 
occur. The new regions of flux are observed to 
form in the vicinity of the footpoints of spine field 
lines in the pre-reconnection field. Together they 
cover a region of comparable scale to the distri¬ 
bution of current and flux rope structures, here of 
order 10 -1 D. 


5.2. Implications for solar wind models 

Our results imply that in the vicinity of open 
spine structures and open separatrix curtain struc¬ 
tures, an efficient mixing of open and closed mag¬ 
netic flux, and the associated plasma, is likely 
to take place whenever reconnection occurs at 
the corresponding nulls or separatrices. This is 
an attractive ingredient for explaining observed 
properties of the slow solar wind by the inter¬ 


change reconnection model. In particular, the 
slow solar wind is known to be highly fluctuating 
in both composition and velocity (in both space 
and time), with the composition properties vary¬ 
ing from close to those of the closed corona to 
nearly photospheric ( Geiss et al.||1995' Zurbuchen 
2007). Contributing factors to this fluctuating, 


filamentary structure could be the bursty nature 
of the interchange reconnection, and the complex 
spatial structuring on large and small scales of 
the open-closed boundary. Combining our results 
with those from simulations of current layer insta¬ 
bilities, it is clear that the reconnection process 
should lead to a highly dynamic magnetic topol¬ 
ogy in which regions of open and closed flux are 
born and evolve in a complex pattern. 

Interchange reconnection models for solar wind 
acceleration share the common feature that they 
require regions of open flux at photospheric 
heights that are at least predominantly sur¬ 
rounded by closed flux. This is consistent with 
observations of significant components of solar 
wind outflow emanating from locations adjacent 


to active regions (e.g. 

Neugebaue 

et al. 

2002 

Woo et al. 2004 

Brooks et al. 2015 

). The models 


of Fisk and co-workers take a statistical approach 
to the evolution of the Sun’s open flux, in which 
they assume a random orientation for the closed 
loops and thus an isotropic diffusion of open field 
lines as they random walk through the closed flux 
region. One piece of observational evidence that 
is argued to support their model of open field dy¬ 
namics is the coincidence of extended open field 
regions with minima in the local rate of flux emer¬ 
gence (|Fisk||2005[ [Abramenko et ak||2006[ [Hage- 


naar et al. 2008). Interestingly, the models also 


predict that the random walk of open field lines 
will induce a braiding of held lines in the helio¬ 
sphere, and this is also a prediction of the MHD 
simulations of Wyper and Pontin (2 014a|b ) (see 
Figure [5]), though in our model the braiding is in¬ 
duced by turbulent dynamics in the reconnecting 
current layer. 

The S-web model, by contrast to the models 
of Fisk and co-workers, seeks to identify explic¬ 
itly the locations of possible interchange reconnec¬ 
tion - and thus outflow - by analysing the detailed 
magnetic field topology. In a given magnetic field 
extrapolation open field channels and patches can 
be identified, however there are indications that 
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the number of disconnected open field regions can 
greatly increase when the resolution of the photo- 


spheric magnetogram is increased (Edwards 20141. 


Such an increase corresponds to an increase in the 
number of arcs in the S-web. In this paper, we 
have argued that when interchange reconnection 
occurs at a structure of the open-closed boundary 

- and the full reconnection dynamics are included 

- it will tend to generate a layer within which open 
and closed flux are mixed. So in this picture, each 
separatrix arc of the S-web becomes not a sin¬ 
gle line but a band within which the open-closed 
flux boundary is highly structured. This may par¬ 
tially mitigate against the concern that the static 
S-web is not space-filling and therefore cannot pro¬ 


vide a continuous slow solar wind outflow (Wang 
et al .|2012 ). The dimensions of this layer of mixed 
open-closed flux are of course crucial, but cannot 
be readily estimated from the present approach. 
They will depend both on the size and geometry 
of the fragmented reconnecting current layer, and 
the overall global field geometry that connects this 
volume out to the heliosphere in any given situa¬ 
tion. While the latter can be estimated from static 
models, elucidating the former will require a full, 
detailed dynamical understanding of the reconnec¬ 
tion process. We note also that our results perhaps 
provide a ‘bridge’ between the Fisk et al. and S- 
web models, in the sense that the broad bands of 
mixed open-closed flux expected to form around 
the arcs of the S-web can be thought of as a regions 
within which the diffusive random walk envisaged 
by Fisk and co-workers becomes highly efficient. 
The indication is then that the broadly uniform 
diffusion coefficient in those models should in fact 
be highly structured in both space and time. How¬ 
ever, to determine the spatial and temporal distri¬ 
bution of these regions of efficient open-closed flux 
mixing - and thus their overall effect will require 
a statistical study of the structures of the open- 
closed flux boundary, together with the results of 
the dynamical simulations suggested above. 

5.3. Morphology and structure of flare rib¬ 
bons 

We now consider the implications of our results 
for signatures of particle acceleration in topolo¬ 
gies involving coronal null points and separators. 
Of course, explaining SEP and flare ribbon obser¬ 
vations requires knowledge of particle behaviour, 


and that remains to be studied. However, the 
above results allow us to make a number of predic¬ 
tions. Consider first impulsive SEP events. A se¬ 
ries of recent observational studies has shown that 
impulsive SEP sources are located in open field re¬ 
gions magnetically well-connected with the target 


(Reames 2013 and references therein). However, 
since their composition is more indicative of closed 
field regions, it has been proposed that they are 
accelerated directly during interchange reconnec¬ 
tion in the low corona (Drake et al. 2009). The 


present results demonstrate that particles accel¬ 
erated (by some means) during the reconnection 
process readily have access to open field lines, since 
all field lines are recursively reconnected from open 
to closed within the envelope of mixed flux de¬ 
scribed above. What’s more, when open and 
closed flux are mixed into sufficiently thin layers, 
the distinction between the two is lost for the par¬ 
ticles - certainly the case if the layers are thin¬ 
ner than the Larmor radius. Furthermore, the dy¬ 
namics in the current layer are most likely to be 
turbulent, and braiding of field lines is expected 
(Wyper and Pontin 2014b), both of which are 


known to lead to enhancement of cross-field par¬ 
ticle transport (e.g. Kontar et al.||20lT and refer¬ 
ences therein). 

Our results, combined with those of [Wyper and| 


Pontin (2014a), also provide insights into the ex¬ 


pected structure of flare ribbons in coronal null 
point and pseudo-streamer topologies. In particu¬ 
lar, it has been shown that the separatrix and QSL 
footprints often map onto the locations of the flare 
ribbons. However, the flare ribbons usually ex¬ 
hibit additional structure, often bright kernel-like 


structures that move along the ribbons (Nishizuka 


et al. 

2009 

Barta et al. 

2011) 


These features 
could correspond to the footpoints of the flux rope 
structures formed during 3D current sheet frag¬ 
mentation, these being associated with bundles of 
efficiently mixed open and closed flux. If this were 
the case, one would expect the motion of the bright 
features to be linked to the velocity of the outflow 
from the reconnection region, multiplied by some 
factor resulting from the geometry of the magnetic 
connection between the reconnection site and the 
photosphere. However, as shown by | Wyper and| 


Pontin (2014b), the flux ropes exhibit a compli¬ 


cated dynamics as they kink and interact with one 
another in the reconnection region, so the motion 
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of their footprint on the photosphere is expected 
to deviate significantly from a simple advection. 
We note that performing the same procedure as 
herein - adding a flux ring to simulate the effect 
of reconnection for a background field defined by 
a QSL or hyperbolic flux tube leads not to more 
structure in the associated Q-map, but simply to 
a break of the high-Q layer (results not presented 
here). This may correspondingly imply that the 
signatures of bursty reconnection in such a field 
topology are different to the case when separatri- 
ces are present. 


One further specific conclusion that can be 
drawn regards the nature of the flare ribbons 
associated with the spine footpoints in the iso¬ 
lated dome topology. Masson et al. (20091 and 


Reid et al. (2012) reported observations of the 


flare ribbons in such a topology and noted that 
the ribbons that were postulated to be connected 
with the spine footpoints were extended structures 
with high aspect ratio (rather than being circu¬ 
lar). They proposed that this was related to the 
distribution of the squashing factor for field lines 
around the spine in the associated extrapolated 
equilibrium field. The results of Section [3] provide 
us with an alternative hypothesis: the extended 
elliptical ribbons may mark out the imprint of the 
vertical separatrix curtains (and surrounding arcs 
of high Q) associated with a fragmented current 
layer. It may also be that both effects are impor¬ 
tant for the spine footpoint ribbon extension. Of 
course it remains to be seen what accelerated par¬ 
ticle distributions are expected during a dynamic 
reconnection event, and this will be pursued in a 
future study. 


6. Summary and conclusions 


Magnetic reconnection in the solar corona is 
likely to occur in highly fragmented current lay¬ 
ers, as demonstrated in recent 3D simulations 


( Daughton et al.|2011 Wyper and Pontin|2014a |. 
Here we have used simple static magnetic field 
models to investigate the implications of the cur¬ 
rent layer fragmentation on the large-scale topol¬ 
ogy of representative solar coronal field structures. 
We have shown that this fragmentation can vastly 
increase the topological complexity beyond that 
of the equilibrium magnetic field. In particular, 
when the fragmenting current layer forms at the 


open-closed magnetic flux boundary, the structure 
of that boundary can become highly complex. The 
results, however, are also relevant for the stud¬ 
ied topologies in the case where all flux is glob¬ 
ally closed - in this case some flux closes at some 
distant point on the photosphere (but is ‘locally 
open’). 

We considered here two principal topologies; 
the isolated null point dome and the separatrix 
curtain topology. Both of these are observed over 
a broad range of characteristic scales in the corona, 
from hundreds of Mm down to tens of km. In 
the presence of an isolated dome, non-linear tear¬ 
ing of the reconnecting current layer leads to the 
formation of an envelope of magnetic flux around 
the initial dome structure in which flux from in¬ 
side and outside the dome is efficiently mixed to¬ 
gether. Magnetic flux is continually recursively 
reconnected from open to closed and back again 
within this envelope. The result for the field 
at large heights is that a flux tube is present 
around the original spine line within which field 
lines are being continually reconnected with those 
from the closed region beneath the dome. Such 
isolated dome structures are typically found in 
abundance in coronal holes in magnetic field ex¬ 
trapolations. In the separatrix curtain (‘pseudo¬ 
streamer’) topology of Section [4] we saw that the 
breakup of the current layer leads to the forma¬ 
tion of new flux domains. In particular, open and 
closed flux (as well as flux from pairs of discon¬ 
nected open field regions) form in nested domains 
with very short length scales. The thickness of 
the adjacent open and closed flux domains can 
be many orders of magnitude smaller than the 
global length scale of the field structure or in¬ 
deed of the flux ropes in the current layer - in 
our models with only three flux ropes the map¬ 
ping layers were two orders of magnitude smaller 
than the flux ropes. The expectation is that in 
a dynamic evolution, continual reconnection be¬ 
tween the narrow layers of open and closed flux 
would occur within a flux envelope surrounding 
the new nested flux domains. Our static models 
predict that in the corona immediately above the 
pseudo-streamer this envelope will cover a region 
of comparable scale to the distribution of current 
and flux rope structures. However, this would be 
expected to widen with height as the field strength 
reduces with radial distance from the sun and the 
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field expands laterally. 

Understanding particle acceleration in topolo¬ 
gies such as those studied could help us compre¬ 
hend both the release of impulsive SEPs to open 
field lines and the appearance of certain flare rib¬ 
bons. In particular, the appearance of extended 
flare ribbons at spine footpoints in the dome topol¬ 
ogy could be related to the separatrix footprints 
that appear there when null point bifurcations oc¬ 
cur as the current layer fragments. What’s more, 
the efficient mixing of open and closed flux in the 
reconnection process provides a natural mecha¬ 
nism for accelerated particles to access the open 
field region. Future studies of particle acceleration 
during the reconnection process will reveal much 
more. 


Global magnetic field extrapolations are now 
revealing the huge complexity of the coronal field, 
and in particular the structure of the boundary be¬ 
tween open and closed magnetic flux. Regions of 
open flux that are either disconnected from the po¬ 
lar coronal holes at the photosphere or connected 
only by narrow open-flux corridors contribute arcs 
to the S-web (Antiochos et al. 2011). Our re¬ 
sults show that whenever reconnection occurs at a 
null point or separator of the open-closed bound¬ 
ary, the associated separatrix arc of the S-web be¬ 
comes not a single line but a band of finite thick¬ 
ness within which the open-closed flux boundary 
is highly structured. The dimensions of this band 
are of course crucial, but cannot be readily esti¬ 
mated from the present approach. The next step 
then, to determine the importance of this effect, 
requires dynamical MHD simulations of the pro¬ 
cess in order to quantify the dimensions of this 
band and the flux associated with it. 
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Table 1: Details of model field configurations. s„ 
is the strength of the n th flux ring ( B 0 in Equation 
[2]), p n is the position of its centre ((xn,Vn, Zn) in 
Equation [ 2 ]) . 
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“Background magnetic field where “dome” corresponds to Equation 0 and “curtain” to Equation ([3]) 
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